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ABSTRACT 

We have acquired 65 h of single-site time-resolved CCD photometry of the pulsat- 
ing DB white dwarf star CBS 114 and 62 h of two-site high-speed CCD photometry of 
another DBV, PG 1456+103. The pulsation spectrum of PG 1456+103 is complicated 
and variable on time scales of about one week and could only partly be deciphered 
with our measurements. The modes of CBS 114 are more stable in time and we were 
able to arrive at a frequency solution somewhat affected by aliasing, but still satis- 
factory, involving seven independent modes and two combination frequencies. These 
frequencies also explain the discovery data of the star, taken 13 years earlier. 

We find a mean period spacing of 37.1 ± 0.7 s significant at the 98% level be- 
tween the independent modes of CBS 114 and argue they are due to nonradial g-modc 
pulsations of spherical degree I = 1. We performed a global search for asteroseismolog- 
ical models of CBS 114 using a genetic algorithm, and we examined the susceptibility 
of the results to the uncertainties of the observational frequency determinations and 
mode identifications (we could not provide m values). The families of possible solu- 
tions are identified correctly even without knowledge of m. Our optimal model suggests 
T cff = 21 000 K and M* = 0.730 M Q as well as log(M He /M») = -6.66, X Q = 0.61. This 
measurement of the central oxygen mass fraction implies a rate for the 12 C (a, 7) 16 O 
nuclear reaction near 6*300 = 180keVb, consistent with laboratory measurements. 

Key words: stars: variables: other - stars: variables: ZZ Ceti - stars: oscillations - 
stars: individual: CBS 114 - stars: individual: PG 1456+103 



1 INTRODUCTION 

The helium-atmosphere DB white dwarf stars are very inter- 
esting from the standpoint of stellar structure and evolution. 
The chemical evolution of their atmospheres cannot be satis- 
factorily explained to date (see Shipman 1997 for a review). 
In particular, the presence of the so-called DB gap, the ab- 
sence of DB white dwarfs between temperatures of ~30000 
to 45000 K (Liebert 1986), is poorly understood. It is also 
not clear whether DBs are mostly produced by single-star 
evolution or whether a significant fraction originate from 
binary progenitors. 

The examination of these problems can be aided by as- 
teroseismology - the study of the interiors of pulsating stars 
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via the analysis of their normal mode spectra. Fortunately, a 
class of pulsating DB white dwarf stars (hereinafter DBVs) 
exists, and their prototype, GD 358, is one of the classical 
examples for the successful application of asteroseismologi- 
cal methods (Winget et al. 1994, Vuille et al. 2000). 

Compared to other classes of pulsating star, the DBVs 
seem quite promising candidates for asteroseismology: their 
mode spectra are shown or believed to be rich, and pulsation 
theory for these stars is quite advanced and well-tested. For 
instance, the central oxygen abundance of the models has 
been shown to have a measurable effect on the pulsation fre- 
quencies (Metcalfe, Winget & Charbonneau 2001), and the 
possible presence of a 3 He/ 4 He transition zone due to chem- 
ical diffusion has also been shown to produce a measurable 
effect (Montgomery, Metcalfe & Winget 2001), although re- 
cent observations have ruled out this possibility (Wolff et 
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al. 2002). In addition, the recently developed application of 
genetic-algorithm-based model-fitting to DBV mode spectra 
(Metcalfe, Nather & Winget 2000) allows an objective and 
effective exploration of parameter space to find an optimal 
model for the observations. 

The application of these methods to DBVs other than 
GD 358 is however hampered by a lack of suitable observa- 
tional data. Most of these stars are rather faint (around 16 th 
magnitude) , which requires the use of 2-m class telescopes to 
acquire suitable data with photomultiplier detectors. How- 
ever, time on these telescopes is generally oversubscribed 
and long runs are difficult to obtain. With the advent of 
CCD detectors and with the development of systems capa- 
ble of acquiring measurements with high time resolution, 
the DBVs are now in reach of smaller telescopes. We have 
therefore started to acquire measurements aiming at making 
as many DBVs as possible accessible to asteroseismology. 

The pulsations of the DB white dwarf star CBS 114 
were discovered by Winget & Claver (1988, 1989), who 
reported multiperiodic oscillations with peak-to-peak light 
variations up to 0.3 mag and with a time scale of about 
650 s. PC 1456+103 (hereinafter referred to as PC 1456) was 
reported as a DBV by Grauer et al. (1988). These authors 
reported multiperiodic pulsations with dominant periodic- 
ities between 420 and 850 seconds and maximum peak-to- 
peak amplitudes around 0.15 mag. As no further observa- 
tions of either of these stars have been published, we deemed 
it worthwhile to obtain larger data sets to study the pulsa- 
tions of CBS 114 and PG 1456 in more detail. 



Table 1. Journal of the observations. The SARA measurements 
were obtained between JD 2451987-2451990; the remainder are 
SAAO data 



CBS 114 PG 1456+103 



Run start 


Length 


N 


Run start 


Length 


N 


JD-2450000 


h 




JD-2450000 


h 




1940.421 


3.32 


403 


1940.570 


1.44 


225 


1941.400 


3.74 


450 


1941.609 


0.55 


87 


1942.394 


4.00 


446 


1942.570 


1.49 


219 


1943.446 


1.30 


143 








1944.427 


3.40 


404 


1944.573 


1.51 


233 


1945.397 


4.25 


484 


1945.577 


1.44 


209 


1946.386 


4.17 


435 


1946.567 


1.75 


267 


1954.379 


0.68 


80 








1955.359 


4.20 


382 


1955.539 


2.50 


410 


1956.357 


4.30 


513 


1956.539 


2.62 


400 


1957.353 


4.25 


495 


1957.533 


2.74 


442 


1959.346 


4.38 


464 


1959.530 


2.93 


449 


1960.346 


4.30 


511 


1960.530 


0.55 


81 








1987.929 


2.71 


251 








1988.898 


3.34 


277 








1989.882 


3.60 


349 


1996.249 


4.17 


497 


1996.426 


5.83 


897 


1997.239 


4.32 


511 


1997.423 


5.88 


909 


1998.235 


4.36 


520 


1998.422 


4.37 


617 


1999.355 


1.37 


156 


1999.416 


4.99 


692 


2001.319 


1.38 


107 


2001.415 


4.42 


681 


2002.252 


2.77 


324 


2002.373 


6.94 


808 


Total 


64.66 


7325 




61.60 


8503 



2 OBSERVATIONS 

Most of our measurements were acquired as differential 
CCD photometry with the 0.75-m telescope at the Suther- 
land station of the South African Astronomical Observatory 
(SAAO), at which we had been allocated three weeks of ob- 
serving time. 

We used the University of Cape Town CCD camera 
(O'Donoghue 1995). Its CCD was operated in full-frame 
mode to acquire a maximum number of local comparison 
stars. We used 2x2 or 3x3 prebinning depending on seeing, 
which results in readout times of about 2-3 seconds. For the 
faint CBS 114 (B ~ 17.2 as estimated from our CCD frames 
compared to PG 1456+103), we used integration times of 27 
or 28 seconds to obtain one frame per 30 seconds, and we 
used 20-second integrations for PG 1456. No filter was used 
to maximize the number of photons detected. On every clear 
night, sky flatfields were taken during twilight. 

In addition, we acquired measurements of PG 1456 with 
the 0.9-m SARA telescope at Kitt Peak National Observa- 
tory (Arizona), to reduce the aliasing problem in the fre- 
quency analysis. The CCD used is an Apogee AP7p, with 
a back- illuminated SITe SIA-502AB 512x512 detector. The 
pixels are 24/im square, or 0.75 arcsec/pixel. Read noise for 
the camera is 12.2 electrons r.m.s., and the gain is 6.1 elec- 
trons per ADU. We used integration times of typically 30 
seconds, and our readout time was about 7 seconds. Again, 
we took sky flats during each twilight. Our measurements 
are summarized in Table 1. 

We started data reduction with the corrections for bias, 
dark counts (SARA measurements only) and flatfield. Mean 



weekly flatfields were used for the SAAO data sets and com- 
bined nightly flats for the SARA frames. Photometric mea- 
surements on these reduced frames were made with the pro- 
gram MOMF (Kjeldsen & Frandsen 1992), which applies a 
combination of PSF and aperture photometry relative to a 
user-specified ensemble of comparison stars. No variability of 
any star other than the targets in the different CCD fields 
was found, and the comparison star ensemble resulting in 
the lowest scatter in the target star light curves was chosen. 

The resulting differential light curves were corrected for 
differential colour extinction and for correlations with see- 
ing. For CBS 114, it was sometimes also found necessary to 
de-trend the data with (x, y) position of the star on the chip 
because of flatfield errors. Residual low-frequency trends in 
the data (not found to be coherent over the time span of the 
observations, hence judged not to be intrinsic to the stars) 
were removed by means of low-order polynomials. 

Finally, the times of measurement were transformed to 
a homogeneous time base. We chose Terrestrial Time (TT) 
as our reference for measurements on the Earth's surface 
and applied a correction to account for the Earth's motion 
around the solar system's barycentre. As this barycentric 
correction varied up to ±1 s throughout a run, we applied it 
point by point. Our final time base therefore is Barycentric 
Julian Ephemeris Date (BJED). The reduced time series 
were subjected to frequency analyses; we show example light 
curves in Figs. 1 and 2. The typical rms scatter per single 
data point is around 40 - 50 mmag for CBS 114, and for 
PG 1456 about 15 mmag (SAAO data) and 11 mmag (SARA 
data) . 
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Figure 1. An example light curve of CBS 114; a multifrcquency 
fit is included as well. Although the data are not of impressive 
quality due to the faintness of the star, the multiperiodic pulsa- 
tions are clearly visible. 



3 FREQUENCY ANALYSIS 

Our frequency analysis was performed with the pro- 
gramme PERIOD 98 (Sperl 1998). This package applies 
single-frequency power spectrum analysis and simultane- 
ous multi-frequency sine-wave fitting, but also has some ad- 
vanced options. In particular, it can be used to calculate op- 
timal solutions for multiperiodic signals including harmonic, 
combination, and equally spaced frequencies, which are of- 
ten found in the analysis of the light curves of pulsating 
white dwarf stars. 



3.1 CBS 114 

A frequency analysis of single-site data of a multiperiodic 
variable in the presence of noise is difficult because of alias- 
ing; extreme caution is required. However, the structure and 
extent of our data set aids us in this effort. 

We started by calculating nightly amplitude spectra of 
our data. The amplitude spectra of the individual nights 
were quite similar, with a number of dominant peaks always 
occurring at similar frequencies. However, some variations 
in the individual amplitudes were noted. 

As the next step, we combined the three weekly data 
sets (which are very similar in length and extent) and we 
computed their amplitude spectra as well as that of all three 
weeks combined. They are shown in Fig. 3, which shows six 
dominant structures between 1.4 and 2.6 mHz and some 
further signals which could be combination frequencies. 

Some trial prewhitening within the strongest features 
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Figure 2. An example light curve of PG 1456+103 with a mul- 
tifrcquency fit. As the star is 1.3 mag brighter than CBS 114, the 
data are of much better quality; the multiperiodic pulsations are 
readily visible as well. 



convinced us that all of them are dominated by one signal 
each. Consequently, we determined their frequencies by ex- 
amining the combined, the weekly and the average of the 
weekly amplitude spectra and by least-squares fitting to the 
light curves. In this way, we could rule out a number of 
alias frequencies because they gave unreasonable results if 
adopted. However, in some cases, no unambiguous decision 
could be made. 

Having determined the best set of frequencies, we re- 
fined their values by calculating an optimal solution for the 
whole data set, prewhitened it from the data, and calculated 
the amplitude spectrum of the residuals (upper panel of Fig. 
4). Because of the temporal changes in the weekly amplitude 
spectra, it is not surprising that some residual mounds of 
amplitude remained near the dominant frequencies. 

We examined the frequencies of these features and 
their aliases and performed trial prewhitening to check for 
multiplet structure possibly present. No such evidence was 
found. We then fit our best frequencies to the individ- 
ual nights longer than 3 h (assuring the dominant struc- 
tures are resolved) and determined their amplitudes, phases, 
and corresponding error estimates (following Montgomery 
& O'Donoghue 1999). We found that significant amplitude 
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Figure 3. Weekly and combined amplitude spectra of our 
CBS 114 data. The same peaks are always present, but with some- 
what variable amplitudes. 

and phase changes occurred over time scales of 3 - 4 days, 
and that either amplitude or phase variability alone is not 
sufficient to explain the observations; both are present. 

As we cannot determine the cause of these variations, 
we treated them phenomenologically and calculated our final 
multifrequency solution with constant frequencies as deter- 
mined before, but with amplitudes and phases variable over 
half-week periods. In this way, we can determine ranges of 
the amplitudes of the dominant signals, and most of the 
residual mounds in the amplitude spectrum are thus re- 
moved (lower panel of Fig. 4). In Table 2, we list the best 
multifrequency solution determined with this method. 

3.1.1 Re- analysis of Winget & Claver's data 

Although the measurements by Winget & Claver (1988, 
1989) only comprise 5.2 h of data obtained during two nights, 
we can still use them for a comparison with our results. The 
amplitude spectrum of these measurements (upper panel of 
Fig. 5) shows dominant features at frequencies agreeing very 
well with the ones determined from our data. Consequently, 
we fitted the first six frequencies from Table 2 to these mea- 
surements and noticed that they were consistent with the 




Frequency (mHz) 



Figure 4. Upper panel: residual amplitude spectrum after 
prewhitening the best multifrequency solution for our measure- 
ments of CBS 114, assuming constant amplitudes and phases over 
the whole data set. Lower panel: residual amplitude spectrum af- 
ter prewhitening the same frequency solution, but allowing for 
amplitudes and phases to vary over half- week intervals. 



Table 2. Multifrequency solution for CBS 114. For frequencies 
labeled with a minus sign, the negative daily alias (/j — 11.60/iHz) 
may be the correct frequency, and for frequencies labeled with a 
plus sign, the positive daily alias (/j + 11.60/iHz) cannot be ruled 
out as the correct value. Amplitudes are listed in the ranges they 
assumed during the observations. 



ID Freq. Period Amp. (2001) Amp. (1988) 



(/iHz) (s) (mmag) (mmag) 



fi 1518.75" 658.43 16-33 33 

h 1613.12 619.91 10-17 15 

h 1729.73 578.13 22-37 18 

h 1969.60 507.71 11-17 < 5 

fs 2306.38+ 433.58 4-12 13 

fa 2509.89- 398.42 4-10 9 

fi + f 3 3248.48 5-7 < 7 

h + h 3342.85 2-8 < 6 

ft 1835.64+- 544.77 < 4.2 13 



older data. We therefore determined corresponding frequen- 
cies, amplitudes and phases, subtracted this fit from the 
data and calculated a residual amplitude spectrum, which 
is shown in the lower panel of Fig. 5, revealing the presence 
of yet another signal. 

A variation at this frequency is also present in our new 
observations, although its amplitude is below our detection 
threshold. Consequently, we include it in our final frequency 
solution for the combined data sets (Table 2) . We note some 
more suspected signals in the lower panel of Fig. 5, but due 
to the small amount of data we do not push the analysis 
further. We have therefore detected altogether seven inde- 
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Figure 5. Upper panel: amplitude spectrum of the measurements 
by Winget & Claver (1988, 1989). Lower panel: residual amplitude 
spectrum after prewhitening the first six frequencies determined 
from our observations. 



pendent pulsation frequencies and two combination signals 
in the available light curves of CBS 114. 



3.2 PG1456 

Due to our two-site coverage of the light variations of 
PG 1456, the frequency analysis should be comparatively 
easy. However, initial trials showed that this star's pulsa- 
tional amplitudes and frequencies are also somewhat vari- 
able, and as our light curves from the first two weeks are 
rather short, they complicate the analysis if the total data 
set is used. Consequently, we used the JD 2451987 - JD 
2452002 two-site measurements as our primary data set for 
the frequency analysis; the remaining data were used for 
comparison purposes and consistency checks. Fig. 6 shows 
some prewhitening steps in these two main subsets of data. 

A signal (/s) not present in the early SAAO observa- 
tions reached considerable amplitude in the combined two- 
site data obtained about 7 weeks later; we note that this 
is not due to an aliasing problem. As a matter of fact, the 
amplitude of this signal was still increasing during the latter 
measurements. It had about 12 mmag in the SARA data, but 
18 mmag in the last set of SAAO measurements, in which 
it also showed a significant phase variation. A combination 
frequency (/a+/s) not present initially also appeared in the 
two-site data. On the other hand, the signal fc was found 
to be constant in amplitude and phase over the whole data 
set, and /a was constant in amplitude but showed a slight 
frequency shift of about 0.4 /iHz between our two subsets of 
data. 

The residual amplitude spectra in the lowest panel of 
Fig. 6 still show a number of interesting peaks, several of 
which may be real. To estimate the richness of the frequency 
spectrum of PG 1456, we performed further prewhitening in 
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Figure 6. Amplitude spectra of our measurements of PG 1456. 
The upper halves of the panels are for the two-site JD 2451987 
- JD 2452002 data, whereas the lower, mirrored, halves originate 
from the preceding single-site observations. Uppermost panel: 
original data. Second panel: data after prewhitening of the dom- 
inant frequency. Another frequency and a combination are de- 
tected only in the two-site data. Third panel: another mode is 
found in all the data. Lowest panel: the residual amplitude spec- 
tra; more signals are clearly present. 
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Figure 7. Prcwhitcncd amplitude spectra of our measurements 
of PG 1456. The arrows in the individual panels indicate further 
possible signals which are prewhitcned in the following panels. 



the two-site data (still checking with the early SAAO mea- 
surements), taking the temporal variations of fs into ac- 
count. The resulting amplitude spectra are shown in Fig. 7. 

It becomes clear that the light variations of PG 1456 arc 
quite complicated and that our data set is not sufficient for 
finding all the frequencies present. A preliminary multifre- 
quency solution is given in Table 3. 

A few comments are still necessary. Signals in the region 
around 1970 /iHz are present in both data sets with similar 
amplitude (e.g., the 1973.4 /iHz signal in Table 3), but their 
frequencies are different by 5 /iHz, even taking aliasing into 
account. The dominant frequency region in the residual sin- 
gle site-data is around 1304 /iHz, but neither this frequency 
nor its aliases are present in the two-site data. It seems that 
only intensive multisite observations will make a good un- 
derstanding of the frequency spectrum of PG 1456 possible. 

The periods we determined are in very good agreement 
with those detected by Grauer et al. (1988). Regrettably, 
a detailed comparison of these periods is impossible as the 
discovery data are no longer available (Grauer, private com- 
munication) . 



4 INTERPRETATION OF THE 
OBSERVATIONAL RESULTS 



Table 3. Preliminary multifrcqucncy solution for PG 1456. 
Whereas we are sure about the frequencies /a to fc, the other 
values may be affected by aliasing. Amplitudes are given for both 
the two-site measurements (left column) and the preceding single- 
site SAAO data. Because of amplitude and frequency variability, 
error estimates are only given for amplitudes of apparently con- 
stant signals. 



ID 


Freq. 


Period 


Amp. 


Amp. 




(MHz) 


(s) 


(mmag) 


(mmag) 




Certain detections 




!a 


1505.9 


664.1 


18.9 ± 0.3 


18.5 ± 0.6 


f B 


1465.4 


682.4 


11.7 - 18.4 




1a + }b 


2971.3 




5.7 - 6.8 




fc 


2387.9 


418.8 


4.0 ± 0.3 


5.0 ± 0.6 




Possible further signals 






1195.9 


836.2 


4.4 ± 0.3 


4.4 ± 0.6 




1973.4 


506.7 


2.7 ± 0.3 






1224.8 


816.5 


3.9 ± 0.3 






1324.9 


754.8 


3.9 ± 0.3 






Period spacing (s) 



Figure 8. Search for a preferred period separation within the 
independent signals in the light curves of CBS 114. A spacing of 
37.1s and its harmonics dominate this diagram. 



4.1 CBS 114 

The independent pulsation periods of CBS 114 as listed in 
Table 2 seem to have similar spacing. We therefore searched 
for and determined their mean spacing by calculating the 
Fourier power spectrum of the period values with unit am- 
plitude (Handler et al. 1997), which we show in Fig. 8. A 
mean period spacing (significant at the 98% level) is indi- 
cated, amounting to 37.1±0.7s. This result is corroborated 
by a Kolmogorov-Smirnov test (not shown) . 

We interpret this preferred period separation as a sign 
of the excitation of a number of pulsational radial overtones 
with the same spherical degree I. Some small deviations 
around its mean value (see Fig. 9 in Sect. 5) are an indi- 
cation of mode trapping. The presence of a preferred pe- 
riod separation is predicted by asymptotic theory (Tassoul 
1980), and was already observed within e.g. the I = 1, m = 
modes of GD 358 (Winget et al. 1994, Vuille et al. 2000), for 
which a mean period spacing of 39.2 s was found. In addi- 
tion, Bradley & Winget (1994) showed that the systematic 
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deviations from this mean spacing were due to mode trap- 
ping. 

We cannot directly identify the £, m values of the modes 
we observed in CBS 114 because of the lack of any rota- 
tional m-mode splitting in our frequency spectra. However, 
it seems very likely that the independent modes of CBS 114 
are all £ = 1: if we saw a mixture of £ values, we would not 
expect to find such a significant mean period spacing. There- 
fore the modes we observed must originate exclusively (or at 
least predominantly) from the same I. They should be £ = 1 
or 2 as the effects of geometrical cancellation (Dziembowski 
1977) are expected to render modes of higher £ photometri- 
cally undetectable. 

The mean period spacing between consecutive overtones 
of a pulsating white dwarf star is a measure of its mass 
(Kawaler 1987). If the modes of CBS 114 were I = 1, their 
mean period spacing would be consistent with the star be- 
ing a bit more massive than CD 358. However, if the modes 
were all £ — 2, CBS 114 would need to have a mass below 
O.3M (see Bradley, Winget & Wood 1993). Such a low 
mass is inconsistent with the spectroscopic gravity of the 
star (Bcauchamp et al. 1999) and would be quite unusual 
considering the mass distribution of the DB white dwarf 
stars (Beauchamp et al. 1996). 

Based on these arguments, we suggest that CBS 114 is a 
DBV star pulsating predominantly in nonradial g-modes of 
spherical degree 1=1. We note that its individual pulsation 
periods, especially the longer period modes, are very similar 
to those of GD 358. 



4.2 PG1456 

The interpretation of the frequency spectrum of PG 1456 is 
more difficult. A search for a mean period spacing such as 
in the previous section did not give a significant result. This 
suggests that the modes we detected are not of the same £ 
or that the star is a fast rotator. The modes of m / would 
then mask the possible patterns of equally spaced periods. 

About half of the indpendent periods we detected or 
suspect are similar to periods of GD 358, but the others are 
not. We are therefore unable to interpret the mode spectrum 
of PG 1456 at this stage. Extensive multisite observations, 
e.g. with the Whole Earth Telescope (Nather et al. 1990), 
are required to understand this star's pulsations. 



5 ASTEROSEISMOLOGY OF CBS 114 

Using the optimization method developed by Metcalfe, 
Nather, & Winget (2000) and Metcalfe, Winget, & Char- 
bonneau (2001), we performed a global search for the op- 
timal model parameters to fit the 7 independent pulsation 
periods of CBS 114 listed in Table 2. The method uses a 
parallel genetic algorithm to minimize the root-mean-squarc 
(rms) differences between the observed and calculated peri- 
ods for models with effective temperatures (T cf f) between 
20,000 and 30,000 K, total stellar masses (M») between 0.45 
and 0.95 Mq, and helium layer masses with — log(MH e /M*) 
between 2.0 and ~7.0. This technique has been shown to find 
the globally optimal set of parameters consistently among 
the many possible combinations in the search space, but it 



requires between ~10 and 4000 times fewer model evalu- 
ations than an exhaustive search of parameter-space to ac- 
complish this (depending on the number of free parameters), 
with a failure rate < 10 -5 . 

We assumed that all of the observed modes had spher- 
ical degree I = 1 (as suggested in Sect. 4.1) and azimuthal 
order m = 0. The latter assumption can bias our determina- 
tion of the particular set of model parameters that produces 
the optimal fit to the data, but if the rotation period is ~1 
day, any set of periods drawn from ra = (— 1, 0, +1) will pro- 
duce essentially the same overall picture of the parameter- 
space. We demonstrated this general behavior by generating 
C-core fits to 100 data sets for GD 358 that used randomly- 
selected m-components from those identified in Winget et 
al. (1994), and the m = values for k = 12, 18 (for which no 
triplet structure was found). In every case, the optimal set 
of model parameters fell within the same families of good 
solutions that were identified using the m = modes. Thus, 
when the spherical degree of the modes is known, the model- 
fitting procedure correctly identifies the families of possible 
solutions even when the values of m are unknown. Further- 
more, we found that the optimal solutions fell into the var- 
ious families within parameter-space in proportion to the 
relative fitness of that family when only m — modes were 
used for the fit. For example, if a family has a peak fitness in 
the m = case that is 3a better than any other family, then 
the use of modes will identify an optimal solution in 

this same family with ~50% probability. This means that if 
one of the families produces much better solutions than the 
others, we can be reasonably confident that we have at least 
identified the correct family of solutions. 

5.1 Carbon-core fit 

Assuming a pure C core extending to a fractional mass of 
0.95 m/M*, the optimal set of model parameters found by 
the genetic algorithm for CBS 114 were 

T off = 24,600 K, M, = 0.655 M , log(M He /M») = -3.96, 

with rms period residuals ap — 1.02 s. The observed and 
calculated periods are shown in the top panel of Fig. 9 plot- 
ted against the deviations from the mean period spacing 
(dP), which were calculated using the same set of periods in 
both cases. Each point in this representation of the data is 
independent of the others, unlike a period spacing diagram 
using AP (= P k+1 - P k ; cf. Bradley & Winget 1994). Note 
that the genetic algorithm only fits the periods of the pulsa- 
tion modes, and the agreement between the deviations from 
the mean period spacing is simply a reflection of the overall 
quality of the match. 

As noted in Table 2, the daily aliases of several of the 
identified modes in CBS 114 cannot be ruled out entirely. 
We can assess the impact of this uncertainty on the final set 
of optimal model parameters by repeating the fitting proce- 
dure using one or more of the alias periods. However, since 
the genetic algorithm fitting method with 3 free parameters 
is only about ten times more efficient than calculating the 
entire grid of models, it is better to calculate the pulsation 
periods of all 10 6 models if we intend to repeat the proce- 
dure more than a few times. As there are four periods with 
a total of five possible aliases, we chose to calculate the en- 
tire grid of models. This allows us to check the answer that 
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Figure 9. The observed periods of CBS 114 (solid) and the opti- 
mal model periods found by the genetic algorithm (open) plotted 
against the deviations from the mean period spacing, assuming a 
core of (a) pure C (b) pure O, and (c) mixed C/O (see text for 
complete details). 
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Figure 10. Front and side views of the search space assuming an 
C-core (top panels) and an O-core (bottom panels). Square points 
mark the location of models that yield a reasonable match to the 
periods observed for CBS 114. The darkness of a point indicates 
the relative quality of the match (see text for complete details). 



resulted from the genetic algorithm fit, and will enable us 
to generate C-core fits very quickly in the future for any 
additional data sets on DBV white dwarf stars. 

All C-core models with rms period residuals smaller 
than 2.22 seconds are shown in the top half of Fig. 10, which 
includes front and side views of the entire search space. Each 
point in the left panel corresponds to a point in the right 
panel, and the darkness of a point indicates the relative qual- 
ity of the match with the observations. Black points are 
within Aap — 0.03 s of the optimal model, and the four 
progressively lighter shades of grey correspond to models 
within 3, 10, 25, and 40 times this difference. The parameter- 
correlations explained by Metcalfe, Nather, & Winget (2000) 
are clearly visible, causing the good models to fall along lines 
in the plot rather than on a single point. 

As expected, the optimal parameters from comparison 
of the observations to the complete grid of models were iden- 
tical to those found using the genetic algorithm method. 
In fact, the identified periods led to lower rms residuals 
than any set of periods that included one of the possible 
aliases. Of the combinations with two periods replaced by 
their aliases, fits with lower residuals were achieved in three 
cases: (1) fi /f and f 7 -> , (2) fi -> /f and fr -> /f , 
and (3) fe — > f$ and fj — > fy. When either three or four 
periods were replaced with aliases, only one combination led 
to lower residuals: /i — ► f±, fs — * f^~, and fr — ► . Note 
that aliases of /i and /V appear in most of these alternative 
period lists, implying that they are relatively important to 
the outcome of the fit. Although lower residuals were pos- 
sible using these various combinations of two or more alias 
periods, the optimal models in every such case are too mas- 
sive and too cool to be reconciled with the spectral line fits 
for CBS 114 by Beauchamp et al. (1999). 



5.2 Oxygen-core fit 

The optimal set of model parameters found by the genetic 
algorithm for a pure O core extending to 0.95 m/M* had 
a mass that was inconsistent with the measurements of 
Beauchamp et al. (1999) and a helium layer mass at the 
edge of our search range, near the theoretical limit where 
nuclear burning will occur at the base of the envelope. The 
second-best model had residuals only ~0.03 seconds higher 
and was similar to the optimal C-core model. The parame- 
ters of this model were 

T cH = 25,800 K, M, = 0.640 M Q , log(M He /M„) = -3.96, 

with rms period residuals op — 0.91 s. The middle panel 
of Fig. 9 shows the calculated periods of this model along 
with the observations, plotted against the deviations from 
the mean period spacing. The fit is only slightly better than 
the C-core model, and is qualitatively similar. 

The O-core models that were calculated by the genetic 
algorithm during the optimization process with rms period 
residuals smaller than 2.08 seconds are shown in the bot- 
tom half of Fig. 10. Although this is not a complete sam- 
pling of the parameter-space, it is heavily sampled in regions 
where models produce better than average residuals. Again, 
the shade of each point indicates the relative quality of the 
match — black points are within Aap = 0.03 s of the optimal 
model, and the shades of grey are at 3, 10, 25, and 40 times 
this difference. Note that because the optimal O-core model 
has residuals ~0.1 seconds lower than the optimal C-core 
model, the darkest grey points in the bottom half of Fig. 10 
are actually better than the the black points in the top half. 
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5.3 C/O-core fit 

Since it may seem dubious to fit the 7 observed periods of 
CBS 114 using a model with five free parameters, we proceed 
cautiously. To allow a systematic exploration of models with 
various internal C/O profiles, Metcalfe, Winget, & Charbon- 
neau (2001) used a simple parameterization that explored 
a general class of profiles similar to those used by Bradley, 
Winget, & Wood (1993). The parameterization fixes the oxy- 
gen mass fraction to its central value (Xo ) out to some frac- 
tional mass (q) where it then decreases linearly in mass to 
zero oxygen at 0.95 m/M*. 

Using this method, the optimal model parameters for 
the observed pulsation periods of CBS 114 were 

T cH = 21, 000 K, M, = 0.730 M , 

log(M Hc /M») = -6.66, Xo = 0.61, q = 0.51, 

with rms period residuals op = 0.43 s. All models within 
Delta sigma = 0.03 s of this fit had the same mass, tem- 
perature, and helium layer mass at our sampling resolu- 
tion of 0.005 Mq, 100 K, and 0.05 dex respectively. Only 
models with a central oxygen mass fraction within Delta 
Xo = ±0.01 and the optimal value of q had rms residuals 
within this range. 

The calculated periods of this model are shown in the 
bottom panel of Fig. 9 along with the observed periods, plot- 
ted against the deviations from the mean period spacing. To 
evaluate whether or not this fit is better by an amount that 
justifies the addition of two free parameters, we can use the 
Bayes Information Criterion (BIC, following Koen & Laney 
2000): 

DTr , j. r [logN ohs \ 2 

where N p is the number of free parameters, iV bs is the 
number of observed periods, and a is the rms period resid- 
ual of the optimal fit. The value of BIC must be lower for 
a decrease in a to be considered significant. Our best 3- 
parameter fit was the O-core model, which leads to a value 
of BIC = 0.28. This leads us to expect the residuals of a 
5-parameter fit to decrease to 0.69 seconds without being 
considered significant. The rms residuals of our 5-parameter 
fit are substantially lower than 0.69 seconds, so the addition 
of the extra parameters seems to be justified, and the C/O 
fit is significantly better than the O-core model. 

Even so, we might worry that our optimal model pa- 
rameters may be less accurate than the values obtained for 
CD 358, due to the smaller number of observed pulsation 
periods. To determine the magnitude of any systematic un- 
certainties due to the smaller number of observed modes, 
we performed a new fit to the pulsation periods of CD 358, 
but using only the 7 modes corresponding to those observed 
in CBS 114: k = 8, 9 and 11-15. When we compared the re- 
sult of this fit with the one using all 11 modes observed in 
CD 358, we found only small shifts to the optimal model 
parameters: 

AT efI = -500 K, AM* = +0.015 M , 

Alog(M Hc /M,) = -0.05, AXo = -0.02, Aq = -0.01. 
By analogy, we might expect our 5-parameter fit to CBS 114 



to represent a slight underestimate of the temperature and 
overestimate of the mass. Both potential biases help to ex- 
plain part of the discrepancy between the mass and tem- 
perature of our optimal model for CBS 114 and the values 
inferred from the spectroscopic analysis of Beauchamp et 
al. (1999) [T cS = 23,300; log g = 7.98]. Additional dif- 
ferences are expected since we have used different mixing- 
length parameters (ML3) than Beauchamp et al. (ML2/a = 
1.25). Metcalfe, Salaris, & Winget (2002) quantified the off- 
sets between fits to CD 358 using ML3 and ML2, and found 
them to be approximately the same size as the shifts due to 
the smaller number of data points. Unfortunately, CBS 114 
does not have a published parallax, so an independent con- 
straint on the mass and temperature from the luminosity is 
not presently available. 



6 CONCLUSIONS 

We have presented time-resolved CCD photometry of the 
pulsating DB white dwarf stars CBS 114 and PG 1456+103. 
Our data, obtained with telescopes of only 0.75 m and 0.9 m 
apertures, are sufficient for the extraction of useful astero- 
seismological information even for a star as faint as CBS 114 
(B « 17.2). We have also shown that it is possible to un- 
derstand the mode spectra of multiperiodic pulsating white 
dwarf stars from single-site observations in suitable cases 
(see Handler 2001 for another example). In other cases, like 
PG 1456, this is not possible, and worldwide observing cam- 
paigns, e.g. with the Whole Earth Telescope, are necessary. 

The frequency analysis of our measurements of CBS 114 
resulted in the discovery of a mean period spacing of 
37.1 ± 0.7 s in the independent modes of the star. Although 
no convincing evidence for rotational splitting within the 
modes was detected, we argued that the star pulsates in 
nonradial I = 1 g-modes. 

We used our data set to check the sensitivity of the 
genetic-algorithm-based model-fitting method (Metcalfe et 
al. 2000) to artifacts of data analysis, like aliasing problems 
and our lack of assignments of the azimuthal order m to 
the observed modes. Encouragingly, our results are hardly 
affected by these shortcomings. 

Most significantly, the genetic algorithm found a C/O 
fit to the pulsation periods of CBS 114 that proved to be 
significantly better than the pure-core fits, with an accu- 
racy comparable to a similar fit to GD 358. The optimal 
mass and central oxygen mass fraction from this fit can pro- 
vide an independent measurement of the astrophysically- 
important rate for the 12 C(a,7) 16 reaction, as was re- 
cently done for GD 358 (Metcalfe, Salaris, & Winget 2002). 
In general, higher mass white dwarfs are expected to have 
a lower central oxygen mass fraction, because the 3a rate 
rises faster with increasing density than the 12 C(a,7) 16 
rate. Our optimal model for CBS 114 has a higher mass 
and a lower central oxygen mass fraction than the optimal 
model for GD 358 (Metcalfe, Winget, & Charbonneau 2001), 
so it is consistent with this general trend. A model of the 
internal chemical profile with the same mass as our fit to 
CBS 114 requires a rate for the 12 C(q,7) 16 reaction near 
5*300 = 180 keV b to produce a central oxygen mass frac- 
tion of 0.61 (M. Salaris, private communication). This value 
is close to the rate derived from recent high-energy labo- 



10 Handler, Metcalfe & Wood 



Wolff, B., Koester, D., Montgomery, M. H., Winget, D. E., 2002, 
A&A, in press 



ACKNOWLEDGEMENTS 

We thank Don Winget and Chuck Claver for permission to 
use their data of CBS 114, Maurizio Salaris for providing a 
new internal chemical profile to match our optimal model 
and Mike Montgomery for carefully reading a draft version 
of this paper. 



REFERENCES 

Bcauchamp, A., Wesemael, F., Bergeron, P., Liebert, J., Saffcr, 

R. A., 1996, in Hydrogen- Deficient Stars, eds. C. S. Jeffery & 

U. Heber, ASP Conf. Scr. 96, p. 295 
Bcauchamp, A., Wesemael, F., Bergeron, P., Fontaine, G., Saffcr, 

R. A., Liebert, J., Brassard, P., 1999, ApJ 516, 887 
Bradley, P. A., Winget, D. E., Wood, M. A., 1993, ApJ 406, 661 
Bradley, P. A., Winget, D. E., 1994, ApJ 430, 850 
Dziembowski, W. A., 1977, Acta Astr. 27, 203 
Graucr, A. D., Bond, H. E., Green, R. F., Liebert, J., 1988, AJ 

95, 879 

Handler, G., et al., 1997, MNRAS 286, 303 
Handler, G., 2001, MNRAS 323, L43 

Kawaler, S. D., 1987, in Stellar Pulsations, eds. A. N. Cox, W. 

M. Sparks & S. G. Starrfield (Berlin: Springer), p. 367 
Kjeldsen, H., Frandsen, S., 1992, PASP 104, 413 
Koen, C. & Lancy, D. 2000, MNRAS, 311, 636 
Kunz, R., Fey, M., Jaeger, M., Mayer, A., Hammer, J. W., Staudt, 

G, Harissopulos, S., Paradellis, T. 2002, ApJ, 567, 643 
Liebert, J., 1986, Proc. IAU Coll. 87, p. 367 
Metcalfe, T. S., Nather, R. E., Winget, D. E., 2000, ApJ 545, 974 
Metcalfe, T. S., Winget, D. E., Charbonncau, P., 2001, ApJ 557, 

1021 

Metcalfe, T. S., Salaris, M., Winget, D. E. 2002 ApJ, in press 
Montgomery, M. H., O'Donoghue, D., 1999, Delta Scuti Star 

Newsletter 13, 28 (University of Vienna) 
Montgomery, M. H., Metcalfe, T. S., Winget, D. E., 2001, ApJ 

548, L53 

Nather, R. E., Winget, D. E., Clemens J. C, Hansen, C. J., Hine, 
B. P., 1990, ApJ 361, 309 

O'Donoghue, D., 1995, Baltic Astronomy 4, 519 

Shipman, H., 1997, in White dwarfs: Proceedings of the 10th Eu- 
ropean Workshop on White Dwarfs, eds. J. Isern et al., Dor- 
drecht: Kluwer Academic Publishers, p. 165 

Sperl, M., 1998, Master's Thesis, University of Vienna 

Tassoul, M., 1980, ApJS 43, 469 

Vuille, F., et al., 2000, MNRAS 314, 689 

Winget, D. E., Claver, C. F., 1988, IAUC 4595 

Winget, D. E., Claver, C. F., 1989, in White dwarfs; Proceedings 
of IAU Colloquium 114, ed. G. Wegner, Springer Verlag, p. 
290 

Winget, D. E. et al., 1994, ApJ 430, 839 



ratory measurements (5*300 = 165 ± 50 keV b; Kunz et al. 
2002). By contrast, the rate derived from the optimal model 
of GD 358 by Metcalfe, Salaris, & Winget (2002) was sig- 
nificantly higher (S300 = 370 ± 40 keV b). This suggests 
either that presently unknown sources of systematic uncer- 
tainty in our models must affect the analysis of GD 358 and 
CBS 114 in different ways, or that the two stars have dif- 
ferent evolutionary origins, or both. An asteroseismological 
determination of the central oxygen mass fraction for ad- 
ditional DBV white dwarfs will help us to decide which of 
these scenarios is most likely. 



